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We study spin glass clusters ("shards") in a random transverse magnetic field, and determine the 
regime where quantum chaos and random matrix level statistics emerge from the integrable limits 
of weak and strong field. Relations with quantum phase transition are also discussed. 
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Quantum manifestations of classical chaos and inte- 
grability have been intensively investigated in the last 
decade Q . It has been realized that the statistical prop- 
erties of the quantum energy spectrum are strongly influ- 
enced by the underlying classical dynamics. Indeed, clas- 
sical integrability generally implies the absence of level 
repulsion and a Poissonian statistics of energy level spac- 
ings. On the contrary, chaotic dynamics leads to level 
repulsion with the level spacing statistics P(s) being the 
same as in the Random Matrix Theory (RMT), i.e. the 
Wigner-Dyson distribution (WD) |||. These two distri- 
butions of P(s) also characterize respectively the local- 
ized and metallic phases in the Anderson model of dis- 
ordered systems. At the critical point between these two 
phases an intermediate level spacing statistics occurs || . 

While for one-particle systems the statistical proper- 
ties of spectra are well understood, the same problem 
in many-body systems has been addressed only recently. 
The first results demonstrated that many-body inte- 
grable systems are still characterized by Poisson statis- 
tics Pp(s), whereas in absence of integrals of motion 
the Wigner-Dyson statistics Pwd(s) has been found |Q. 
More recently, investigations of finite Fermi systems such 
as the Ce atom || and the 28 Si nucleus || put forward 
the question of the statistical description and thermal- 
ization induced by interaction. A quantum chaos crite- 
rion for emergence of RMT statistics and dynamical ther- 
malization induced by interaction was established in . 
Namely, a crossover from Poisson statistics to the WD 
distribution takes place when the coupling matrix ele- 
ments become comparable to the level spacing between 
directly coupled states. For two-body interaction, the 
critical interaction strength becomes exponentially larger 
than the multiparticle level spacing. Indeed the latter de- 
creases exponentially with the number of particles, while 
the former decreases typically only quadratically. This 
result was corroborated in [§) , and should apply to vari- 
ous physical systems such as complex nuclei, atoms, clus- 
ters and quantum dots. 

In this Letter we develop and apply these concepts 
to disordered spin systems. These systems are of great 
experimental and theoretical interest 0]; in particular, 
quantum spin glasses have recently attracted a great deal 
of attention and various approaches have been devel- 



oped to investigate their properties J1Q| . Such important 
problems as zero temperature quantum phase transition, 
structure of the phase diagram, correlations and suscepti- 
bility properties have been intensively investigated both 
analytically and numerically |]ll|-^6|. However, to the 
best of our knowledge the level spacing statistics has not 
been studied in disordered spin systems. Here we present 
theoretical estimates for spin glass " shards" which deter- 
mine the crossover from integrability to quantum chaos in 
a way analogous to the case of finite interacting fermionic 
systems. Numerical investigations of P(s) statistics give 
a powerful test of this crossover both at high tempera- 
ture T and at T — where a quantum phase transition 
is believed to take place. 

The spin glass shards we study are described by the 
Hamiltonian: 
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where the <7j are the Pauli matrices for the spin i and 
the first sum runs over all spin pairs. The local random 
magnetic field is represented by T,; uniformly distributed 
in the interval [0, T]. The exchange interactions Jjj are 
distributed in the same way in [—J/y/n, J/y/n\ where n 
is the total number of spins. For r = this system is 
the classical Sherrington-Kirkpatrick spin glass model M . 
The y/n factor in the definition of Jij ensures a proper 
thermodynamic limit at n — > oo. In the limit J/T — > 
one obtains a paramagnetic phase with all spins in the 
field direction. The opposite limit (J/T — > oo) corre- 
sponds to a spin glass phase at low temperature || . 

Let us first discuss the properties of the model for 
highly excited states near the center of the energy band. 
At J = the system is integrable, since there are n in- 
tegrals of motion, and P(s) should follow the Poisson 
distribution. When J/T increases, we expect integra- 
bility to be destroyed and therefore a crossover towards 
Pwd(s) statistics typical of RMT should take place. At 
r = the model again becomes integrable since there 
are n operators (erf) commuting with the Hamiltonian, 
hence P(s) = Pp(s). The latter may seem surprising, 
especially in the light of recent discussions on chaos in 
classical spin glass jL7|. However, in spite of a possible 
complex thermodynamic behavior (Monte-Carlo dynam- 



1 



ics) , the Hamiltonian dynamics of a classical spin glass is 
integrable. As a result, we expect another crossover from 
a WD statistics back to a Poissonian one for large J/T. 
This picture is confirmed by the numerical results for 
P(s) displayed in Fig.l. These P(s) distributions were 
obtained for the states of the same symmetry class (5*2), 
namely the number of spins up is always an even number 
(the interaction does not mix states with even and odd 
number of spins up). We will call the other symmetry 
class with an odd number of spins up S\. 




FIG. 1. Crossover from Poisson to WD statistics in the 
model (^]) for the states in the middle of the energy band 
(±12.5% around the center) for n=12: J = 0,r] = 0.984 (+); 
J/P = Jcp/F = 0.38,r? = 0.3 (x); J/r = 0.866,77 = 0.027(*); 
J/r = J cs /r = 6.15,7? = 0.3 (o); r = 0,7? = 0.99 (□). Full 
curves show the Poisson and WD distributions. Total statis- 
tics (NS) is more than 3 x 10 4 ; s is in units of mean level 
spacing. Insert shows P(s) for the first excitation from the 
ground state in the chaotic regime for n = 15, J/r = 0.465, 
7? = 0.018 (NS = 3000) (x); the full line shows Pwd(s). 

To analyze the evolution of the P(s) distributions 
with respect to the coupling J, it is convenient to use 
the parameter r? = f*°(P(s) - P WD (s))ds/ Jq"(Pp(s) - 
P\VD{s))ds, where s — 0.4729... is the intersection point 
of Pp(s) and P\vd(s) @- In this way t? = 1 corresponds 
to the Poissonian case, and 77=0 to the WD distribution. 
The data of Fig.l show that the distributions with the 
same value of 7? are very close, even if J/T varies more 
than ten times. It is convenient to determine a critical 
coupling strength J c at which the crossover from Pp(s) 
to Pwd{s) takes place by the condition r](J c /T) = 0.3. 
The variation of 77 with respect to J is shown in Fig. 2 
for different n. The global behavior of 77 is in agreement 
with the above picture of transition between integrabil- 
ity and quantum chaos. Indeed, when J/T increases, 7? 
drops to zero and then starts to grow again back to one 
when the spin glass term in ([[]) dominates. The fact that 
a WD statistics sets in is not trivial. Indeed, in dimen- 



sion d — 1 (the model (|l]) with nearest-neighbor cou- 
pling only and Ju+\ drawn randomly from [0, J], which 
was studied in JT^JT^] ) the statistics remains Poissonian 
(7/ = 1) at arbitrary J/r, as is shown in Fig. 2. The 
critical point J/T = 1 |ll|] does not manifest itself. The 
physical reason for this behavior is based on the fact that 
this d — 1 model can be mapped into a model of non- 
interacting fermions (see e. g. pH). Therefore the total 
energy is the sum of one-particle energies that generically 
leads to Pp(s). As a result RMT can never be applied 
to this model and dynamical (interaction induced) ther- 
malization never takes place. Thermalization can only 
appear through a coupling to a thermal bath that was 
implicitly used in ||l8|.16|. On the contrary in the sys- 
tem ([!]) there are two transitions from integrability to 
chaos which determine two critical couplings J cp , from 
the paramagnetic side, and J cs from the spin glass side 
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FIG. 2. Dependence of 7? on the rescaled coupling strength 
J/ J cp for the states in the middle of the energy band for 
n = 7( V ), 8(D), 9(o), 10(o), ll(+),12(*),13(x), 14(A) and 
symmetry S 2 ; NS varies between 12500 and 160000. The 
full line for n = 12 shows the global behavior. Data for the 
d = 1 model (see text) are given for n = 12 (full circles) as 
a function of J/T (upper scale); NS = 30000. The insert 
magnifies the region near J/J cp = 1. 

In analogy with finite interacting fermionic systems R] 
we expect that quantum chaos sets in when the coupling 
strength U becomes comparable to the spacing between 
directly coupled states A c (U ~ A c ). For small J in the 
middle of the spectrum one has U ~ J/y/n and A c ~ 
16r/n 2 since each state is coupled to n(n — l)/2 states 
in an energy band of ST. This gives the quantum chaos 
border from the paramagnetic side: 



J c 



C p T/n 



3/2 



(2) 



where C p is some numerical constant. 



We note that J cp is 



exponentially larger than the multiparticle spacing A n ~ 
2ttT/2™. For large J the transitions between unperturbed 
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states at J 3> T are determined by the T term in (|j) , and 
have typical value T. The number of such transitions, in 
a typical energy interval J, is n, hence A c ~ J/n. This 
gives the quantum chaos border from the spin glass side: 

J cs « C s Tn, (3) 

where again C s is some constant. 

These theoretical predictions are confirmed by the nu- 
merical results presented in Figs. 2-3 which give C p ~ 16 
and C s « 0.5. We attribute the deviation from (|J) in 
Fig. 3 seen for small n to the fact that for these values 
J cp is rather large; this can slightly modify the unper- 
turbed spectrum and the estimates for A c . Also for small 
n the proximity of the second transition at J cs can affect 
the actual value of J cp . The data shown in the insert of 
Fig. 2 indicate that 77 depends only on the ratio J/J cp ; the 
global scaling behavior is less visible than in the fermionic 
model It]], apparently due to the reasons above. 
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FIG. 3. Critical coupling strength J c as a function of n: 
J cp at band center for symmetry S2 (full circles), and near 
the ground state (full squares: Si; full diamonds: S2 ); Jcs at 
band center for S2 (o), and near the ground state (□: Si; O: 
S3). Full lines show the theory (Eqs. (|), (|)) for the band 
center; the dashed line indicates the quantum critical point 
at T — 0, J/r = 0.5. Logarithms are decimal. 

The analysis above at the band center corresponds to 
a very high energy or temperature. Therefore the prop- 
erties of ([!]) near the ground state energy E g should be 
studied separately. To do that we investigated P(s) for 
the first excitations from E g within the same symmetry 
class. For the class Si, the distribution P(s) was com- 
puted for the first four level spacings starting from E g by 
averaging over different disorder realizations. The first 
excitation and its 77 value (770) was analyzed separately 
while the other three were used to generate one cumu- 
lative P(s) with an 7/ noted by rji ]l9| . There are sev- 
eral physical reasons for this separation: in the presence 
of a gap 770 is related to the gap fluctuations while 771 



characterizes the quasi-particle excitations; in a metallic 
quantum dot with non-interacting electrons 770 reflects 
the chaotic character of the dot giving 770 = while the 
spectrum of higher excitations becomes close to Pp(s) 
with rji w 1 due to the absence of interaction. 




FIG. 4. Dependence of r? on J = J/ (J + T): n = 7 (full 
circles), n — 9(o),n = ll(+),n = 13(x),n = 15 (full trian- 
gles), n = 17(o), n = 19(Q); 2000 < NS < 30000. Full curves 
connect data for n = 7, 17. Insert shows 771 in the region near 
J/r = J q /T » 0.5 in more detail; 6000 < NS < 90000. Data 
are for Si symmetry. 

The global behavior of 770 and 771 as a function of 
J = J/(J + r) is shown in Fig. 4. This parametriza- 
tion naturally represents the two integrable limits J = 
and r = in a symmetric way. In both extreme cases 
(J = 0;1), 770 and 771 approach the Poissonian value 1. 
In between there is a pronounced minimum with 770 ~ 
and 771 w 0.1. When the number of spins in the shard in- 
creases, the values of 77 to the right from the minimum go 
up markedly, whereas the minimal 77 value itself changes 
only weakly. These 77 data suggest the existence of a 
crossing point around J = J q m 0.5T. For J < J q the 
difference between curves for different n is small, whereas 
it is much bigger for J > J q . The data in Fig. 3 show the 
variation of J cp , J cs (for which 771 = 0.3) with n; the be- 
havior of J cs is consistent with J cs — > J q for n — > 00 
for both symmetry classes. Also according to these data 
J cp oc n~ a with < a < 0.2. This value of a is smaller 
than the one expected from the quantum chaos criterion, 
used to derive Indeed, for small J near the ground 
state A c ~ r/n whereas the coupling is still J j^fn, and 
therefore one expects a — 0.5. For large n this would 
imply that 77 drops quickly to zero with J and then in- 
creases again after the crossing point J q . In such a sce- 
nario, in the thermodynamic limit (n — > 00) we expect 
Jcp —> and J cs — * J q , so that 77 changes from zero 
(J < Jq) to one (J > J q ) with some intermediate statis- 
tics at the critical point J q . The determination of the 
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rj value at that point for n — > oo requires larger system 
sizes. We expect that this critical point corresponds to 
the quantum phase transition between paramagnetic and 
spin glass phases discussed in |l(],[l^,[l5|] for a case when 
all I\ = r. The above scenario is similar to the situation 
for the Anderson transition in d = 3 S: as for the An- 
derson insulator, the lack of space ergodicity in the spin 
glass phase implies 77 = 1. However, in the paramagnetic 
phase the situation may be more complicated. Indeed, 
the criterion U ~ A c which gives a — 0.5 assumes that 
the couplings remain small and do not strongly modify 
the unperturbed A c . But near E g the typical mean field 
acting on a spin is of the order J and for J m T/y/n is 
much larger than the local fields I\ ~ T/n. The strong 
mean field near E g can change the effective A c . Such a 
fact was seen for the Ce atom ||. This may lead to an- 
other scenario in which in the thermodynamic limit J cp 
tends to some non-zero value smaller than J q . Such a 
behavior would be in agreement with the small variation 
of J cp with n in Fig. 3. 
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FIG. 5. Dependence of 77 on 5E/T for n = 12 (1000 realiza- 
tions of disorder): J/T = 0.52(*)(lower scale), J/T = 3.46(+) 
(upper scale). Data are for Si symmetry. 

The discussions above dealt with the two extreme cases 
of energy values. The variation of 77 between these two 
limits is shown in Fig. 5 near the critical point and in the 
spin glass phase as a function of the excitation energy 
SE = E — Eg. Both cases show an unusual behavior when 
77 initially grows with energy and starts to decrease only 
later. This tendency is more pronounced near the critical 
point. A possible reason for this behavior is that the rel- 
ative influence of the mean field becomes weaker as 8E/T 
increases. It is also possible that the situation is similar 
to the case of a chaotic quantum dot with noninteract- 
ing electrons discussed above: the ground state is chaotic 
but the interaction between the first quasi-particle exci- 
tations is weak and so initially 77 grows with SE. Larger 
system sizes are required to clarify this issue. 

In conclusion, our study shows an unusual and interest- 



ing interplay between integrability, quantum chaos and 
phase transitions in disordered spin systems. 
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